Quantum reconstruction of an intense polarization squeezed optical state 
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We perform a reconstruction of the polarization sector of the density matrix of an intense polariza- 
tion squeezed beam starting from a complete set of Stokes measurements. By using an appropriate 
quasidistribution, we map this onto the Poincare space providing a full quantum mechanical char- 
acterization of the measured polarization state. 
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Efficient methods of quantum-state reconstruction are 
of the greatest relevance for quantum optics. Indeed, 
they are invaluable for verifying and retrieving informa- 
tion. Since the first theoretical proposals [1| and the 
pioneer experiments determining the quantum state of 
a light field Q, this discipline has witnessed significant 
growth Laboratory demonstrations of state tomog- 
raphy are numerous and span a broad range of physi- 
cal systems, including molecules ions atoms 0], 
spins Q, and entangled photon pairs [||. 

Any reliable quantum tomographical scheme requires 
three key ingredients |9| : the availability of a tomograph- 
ically complete measurement, a suitable representation of 
the quantum state, and a robust algorithm to invert the 
experimental data. Whenever these conditions are not 
met, the reconstruction becomes difficult, if not impossi- 
ble. This is the case for the polarization of light, despite 
the fact that many recent experiments in quantum optics 
have been performed using polarization states. The ori- 
gin of these problems can be traced back to the fact that 
the characterization of the polarization state in terms of 
the total density operator is superfluous because it con- 
tains not only polarization information. This redundancy 
can be easily handled for low number of photons, but be- 
comes a significant hurdle for highly excited states. An 
adequate solution has been found only recently: it suf- 
fices to reconstruct only of a subset of the density ma- 
trix. This subset has been termed the "polarization sec- 
tor" [l(| (or the polarization density operator [ll|) since 
its knowledge allows for a complete characterization of 
the polarization state flil ]. 

The purpose of this Letter is to report on the first theo- 
retical and experimental reconstruction of intense polar- 
ization states. Specifically, we focus on the case of intense 
squeezed states to confirm how, even in this bright limit, 
they still preserve fingerprints of very strong nonclassical 
behavior. 

We begin by briefly recalling some background mate- 
rial. We assume a two-mode field that is fully described 



by two complex amplitude operators, denoted by &h 
and ay, where the subscripts H and V indicate horizon- 
tally and vertically polarized modes, respectively. The 
commutation relations of these operators are standard: 
[oj, a k ] = djk, with j, k 6 {H, V}. The description of the 
polarization structure is greatly simplified if we use the 
Schwinger representation 

•h = \{,a} H a v + Oya H ) , h = %(a H a' v - cJ H a v ) , 



J3 = \(p} H CLH 



v a v) 



(1) 



together with the total photon number N 
ayav- These operators coincide, up to a factor 1/2, with 
the Stokes operators, whose average values are precisely 
the classical Stokes parameters. One immediately finds 
that J = (Ji, J2, J3) satisfies the commutation relations 
distinctive of the su(2) algebra: [J\, J2] — iJ3> and cyclic 
permutations. This noncommutability precludes the si- 
multaneous precise measurement of the physical quanti- 
ties they represent. 

The Hilbert space TL describing the polarization struc- 
ture of these fields has a convenient orthonormal basis in 
the form of the Fock states for both polarization modes, 
namely \nn, nv)- However, it is advantageous to use the 
basis I J, to) of common eigenstates of J 2 and J3. Since 
J = N/2, this can be accomplished just by relabeling 
the Fock basis as \ J, m) = \uh = J + ni, ny = J — m). 
Here, for fixed J (i.e., fixed N), m runs from — J to J 
and these states span a (2J -I- l)-dimensional subspace 
wherein J acts in the standard way. Since any polariza- 
tion observable has a block-diagonal form in this basis, 
it seems appropriate to define the polarization density 
operator as 



Q : 



j 



Qn 



\J,m){J,m'\. (2) 



J— m,m' — — J 



For a reconstruction of the quantum state one first has 
to extract the required data from the tomographical mea- 
surements. The overall scheme of our experimental setup 
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probabilities 



wi(n) = Tr[f5II^(n)] - n (J, m\g\ J, m) n , 



(5) 



where fi^(n) = 7l(n) fi„ 7^(n) and | J, m) n is the eigen- 
state of n • J relative to the eigenvalue m, which is pre- 
cisely a SU(2) coherent state [14j. The final tomogram 
reads 



FIG. 1: (Color online) Setup for efficient polarization squeez- 
ing generation and the corresponding Stokes measurement ap- 
paratus. 



is illustrated in Fig. 1. The field to be characterized is 
analyzed using a general polarization measurement ap- 
paratus consisting of a half-wave plate (A/2,8) followed 
by a quarter-wave plate (A/4, (f>) and a polarizing beam 
splitter (PBS). The wave plates transform the input po- 
larization allowing the measurement of different Stokes 
parameters by the projection onto the PBS basis \J,m). 
The PBS outputs are measured directly using detectors 
with custom-made InGaAs photodiodes (98% quantum 
efficiency at DC) and a low-pass filter (< 40 MHz) to 
avoid AC saturation due to the laser repetition rate. The 
RF currents of the photodetectors were mixed with an 
electronic local oscillator at 17.5 MHz and digitized with 
an analog/digital converter at I0 7 samples per second 
with a 16 bit resolution and 10 times oversampling. The 
quantum state we measure is defined by the resolution 
bandwidth of 1 MHz at the 17.5 MHz sideband relative 
to the 200 THz carrier. Ten digitized sample corresponds 
to the photocurrent at this sideband generated by pho- 
tons impinging on the photodiode for 1 /is. This photode- 
tection can be modelled by the positive operator-valued 
measure (POVM) FjJ] 



n; 



I J, m)(J, m\ 



(3) 



so that = Tr(gfl^) is the probability of detecting 
tih = J + m photons in the horizontal mode and simul- 
taneously ny — J — m photons in the vertical one. When 
the total number of photons 2 J is not measured and only 
the difference 2m is observed, the POVM is 



n m = E \J> m )(J> m \ 

J=\m\ 



(4) 



The wave plates in the measurement perform linear 
polarization transformations. These can be described 
in terms of J2, which generates rotations about the di- 
rection of propagation, and J3, which generates phase 
shifts between the modes. In other words, their ac- 
tion is represented by 7^(n) = e l8j2 e^ 3 , where n = 
(cos ^sin#, sin </>sin#, cos 9) is a unit vector given by the 
spherical angles (0,(f>). The experimental histograms 
recorded for each n then correspond to the tomographic 



t (n) = Tr[g n m (n)] = 



E 

J—\m 



a (J,m\g\J,m) n . (6) 



The reconstruction in each (2J + l)-dimensional in- 
variant subspace can be now carried out exactly since it 
is essentially equivalent to a spin J [15]. In fact, after 
some calculations one finds that the tomograms can be 
represented in the following compact form 



w ™( n ) 



1 

2^ 



duj Tr ( g j e 



(7) 



which is precisely the Fourier transform of the character- 
istic function of the observable n • J. 
Inverting this expression, one obtains 



QJ 



2J+1 
dn e 

52 



2tt 

did sin' 
J 

m— — J 



wi(n) e v 



(8) 



where the integration over the solid angle dn = sin 6 < 
extends over the unit sphere 1S2 ■ By summing over all the 
invariant subspaces J, the density matrix can be recon- 
structed as follows 



E 



m— — 00 



dn w m (n) K,(m — n ■ J) 



where the kernel lC(x) is 



lC(x) = 



2J+1 
4tt 2 



2tt 



dto sin ( — 



(9) 



(10) 



We stress the appealing analogy of Eq. © with the more 
widely known formula for the reconstruction of the den- 
sity matrix of a single-mode radiation field from the ho- 
modyne tomograms of the rotated quadratures [H . 

From the exact solution §2$), one can calculate any po- 
larization quasidistribution jl 61 ] - From a computational 
point of view reconstructing the SU(2) Q function turns 
out to be the simplest, since in each invariant subspace 
it reduces to 



Q{J, n) = n (J, m\g\ J, m) n , 



(11) 



so, in view of the form ([5]), it is especially suited for 
our purposes. The evaluation of the Wigner function can 
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also be carried out, although with additional effort. Nev- 
ertheless, we do not expect these two quasidistributions 
to differ notably for the states we study here. As a conse- 
quence, we only need to evaluate the matrix elements of 
the kernel K.(m — n' J). This can be accomplished using 
several equivalent techniques, although the most direct 
way to proceed is to note that 



a {J, m\JC(m — n' • J)| J, m) 



2J+ 1 



4tt 2 

(?H sin (?) cos * 



dtv sin 

2 J 



(12) 



where cosx = n n'. In the limit of J ^ 1 the integral in 
(fT2")l reduces to d 2 8(x) jdx 1 evaluated at x = m — J n • n'. 
Since m can be taken as a quasicontinuous variable, we 
can integrate by parts 

Q(J,n) = ^±l p dm J rfn' rf2 ^ n) S(m- J n-n') . 

(13) 

Thus, in the limit of high photon numbers the reconstruc- 
tion turns out to be equivalent to an inverse 3D Radon 
transform fl7j | of the measured tomograms, which greatly 
simplifies the numerical evaluation of Q(J, n). 

To test this theory, we prepared polarization squeezed 
states by exploiting the Kerr nonlinearity experienced by 
ultrashort laser pulses in optical fibers [13]. Our exper- 
imental setup shown in Fig. 1 uses a Cr 4+ :YAG laser 
emitting 140 fs FWHM pulses at 1497 nm with a rep- 
etition rate of 163 MHz. Using the two polarization 
axes of a 13.2 m birefringent fiber (3M FS-PM-7811, 
5.6 /jm mode-field diameter), two quadrature squeezed 
states are simultaneously and independently generated. 
The emerging pulses' intensities are set to be identical 
and they are aligned to temporally overlap. The fiber's 
polarization axes exhibit a strong birefringence (beat 
length 1.67 mm) that must be compensated. To mini- 
mize post-fiber losses, we precompensate the pulses in an 
unbalanced Michelson-like interferometer that introduces 
a tunable delay between the polarizations [lj|. A small 
part (0.1 %) of the fiber output serves as the input to a 
control loop to maintain a ir/2 relative phase between the 
exiting pulses, producing a circularly polarized beam. 

Since the Kerr effect is photon-number conserving, the 
amplitude fluctuations of the two individual modes H 
and V are at the shot-noise level. This was checked using 
a coherent beam from the laser and employing balanced 
detection. The average output power from the fiber was 
13 mW which, with the bandwidth definition of our quan- 
tum state, corresponds to an average number of photons 
of 10 11 per 1 /is. The Kerr effect in fused silica gener- 
ates squeezing up to some Terahertz, making the choice 
of the sideband in principle arbitrary. The wave plates 
in the measurement apparatus were rotated by motorized 
stages. These scanned one quarter of the Poincare sphere 
in 64 steps for 4> (0°- 22.5°) and 65 steps for 6 (0°-45°), 




FIG. 2: (Color online) Three-dimensional isocontour sur- 
face plot of the measured Q(J, n) function for a polarization 
squeezed state. In the inset, we show the projections over the 
coordinate planes passing through the origin of the ellipsoid. 
The projection on the J\- J3 plane shows typical artifacts from 
the Radon transform. 



a measurement of which took over 5 hours. The rest of 
the data can be deduced from symmetry. For each pair 
of angles, the photocurrent noise of both detectors af- 
ter the PBS was simultaneously sampled 2.9 x 10 6 times. 
Noise statistics of the difference of the two detectors' pho- 
tocurrents were acquired in histograms with 2048 bins, 
resulting in the tomograms w m (n) . The optical intensity 
incident on both detectors was recorded as well. 

To reconstruct the Q function we performed a 3D in- 
verse Radon transform. In Fig. 2 we show the result of 
the reconstruction for a polarization squeezed state. Here 
an isocontour surface of the 3D space (or Poincare space) 
of Q(J, n) is seen, as well as the projections of this sur- 
face. The ellipsoidal shape of the polarization squeezed 
state is clearly visible. The antisqueezed direction of the 
ellipsoid is dominated by excess noise stemming largely 
from Guided Acoustic Brillouin Scattering (GAWBS), 
which is characteristic for squeezed states generated in 
optical fibers. Note that this reconstruction of an in- 
tense Kerr squeezed state is very hard with conventional 
homodyne detection techniques for quadrature variables 
due to the high intensities. The projection of the ellip- 
soid on the plane J2-J3 results in an ellipse (and not 
a circle) because of imperfect polarization contrast in 
the measurement setup. As the classical excitation of 
the state is in the J2 direction, one expects to reach the 
shot-noise limit in this projection. After the interference 
at the PBS all intensity impinges on only one detector. 
However, with limited polarization contrast some parts 
of the > 25 dB excess noise from the dark output port 
can be mixed into the intense beam. Acting as a homo- 
dyne measurement this will be visible in the noise of the 
detector. 

In Fig. 3 we compare the isocontour surface plots of a 
coherent and a polarization squeezed state for the value 
corresponding to half width at half maximum. The size 
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FIG. 3: (Color online) Meridional sections of the isocontour 
surface plots of the Q function for a coherent state (blue) and 
a polarization squeezed state (red). 
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of the contours agree with the 6.2 ± 0.3 dB squeezing 
that was directly measured with a spectrum analyzer. 

By summing over all the values of J, we can obtain 
the total Q(n), which is a probability distribution over 
the Poincare unit sphere and is properly normalized. In 
Fig. 4 we have plotted such a function for the squeezed 
state. As the state has a large excitation of the order of 
10 11 photons and the angles of the distribution on the 
unit sphere are small, the spherical coordinates can be 
treated like Cartesian coordinates in the vicinity of the 
classical point and we present a zoomed version of the 
surface of the sphere. Again the excess phase noise of 
the squeezed state is visible. The oscillations with nega- 
tivities near the main peak are characteristic of artifacts 
arising from the inverse Radon transformation. These 
artifacts are more visible for sharp or elongated struc- 
tures. To improve the reconstruction one would have to 
go to more sophisticated reconstructions, e.g., maximum 
likelihood methods 




FIG. 4: (Color online) Probability distribution over the 
Poincare unit sphere for a polarization squeezed state ob- 
tained over all invariant subspaces. The distribution is 
strongly concentrated at the classical mean value, so we show 
a zoomed version. 

In summary, we have presented an exact inversion for- 
mula for quantum polarization states and derived a sim- 
plified version for high-intensity states. The reconstruc- 
tion of an intense polarization squeezed state, formed by 
the Kerr effect in an optical fiber, was performed. Inter- 
esting future investigations include the comparison with 
the maximum likelihood method and the reconstruction 
of nonclassical polarization states with lower intensity. 
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